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• Premise of the study: Vigna unguiculata is an economically important legume, and the complexity of its variability and evolu- 
tion needs to be further understood. Based on publicly available databases, we developed chloroplast microsatellite primers to 
investigate genetic diversity within V. unguiculata and its related species Phaseolus vulgaris. 

• Methods and Results: Twelve polymorphic chloroplast microsatellite markers were developed and characterized in 62 V. un- 
guiculata individuals. The number of alleles per locus varied between two and four, the unbiased haploid diversity per locus 
ranged from 0.123 to 0.497, and the polymorphism information content varied from 0.1 14 to 0.369. In cross-species amplifica- 
tions, nine of these markers showed polymorphism in 29 P. vulgaris individuals. 

• Conclusions: The newly developed chloroplast microsatellite markers exhibit variation in V. unguiculata as well as their trans- 
ferability in P. vulgaris. These markers can be used to investigate genetic diversity and evolution in V. unguiculata and 
P. vulgaris. 
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Cowpea (Vigna unguiculata (L.) Walp.) (2n = 2x = 22), a le- 
gume crop of economic importance, is widely distributed in the 
arid and semiarid regions of Africa, Asia, Europe, Latin America, 
and some parts of the United States (Citadin et al., 2011). As a 
member of the legume family, it belongs to Phaseoleae, the same 
tribe as common bean (Phaseolus vulgaris L.). Compared to its 
close relatives and many other crop species, V unguiculata 
shows a greater tolerance to drought and has the ability to fix ni- 
trogen in poor soils (Muchero et al, 2009). Its grains are a major 
source of dietary protein for humans, and cowpea hay is fed to 
livestock as a nutritious fodder (Badiane et al., 2012). How- 
ever, even though restriction fragment length polymorphism 
(RFLP), amplified fragment length polymorphism (AFLP), 
random amplified polymorphic DNA (RAPD), and simple se- 
quence repeat (SSR) molecular makers have been developed for 
the cowpea nuclear genome, knowledge of variability and evolu- 
tion in the chloroplast genome of V. unguiculata is limited at the 
molecular level (Pro van et al, 2001; Xu et al, 2010). 

Chloroplast microsatellite, or chloroplast simple sequence 
repeat (cpSSR), markers can be used to detect DNA variability 
in the chloroplast genome. They have the same characteristics 
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as nuclear microsatellites, including a multiallelic and codomi- 
nant nature. Moreover, cpSSR markers are found to be poly- 
morphic and transferable among related species because the 
flanking regions of cpSSR loci are conserved. Of particular im- 
portance, cpSSR markers are maternally inherited in most an- 
gio sperms, which allow monitoring of influence on population 
structure by seed-mediated gene flow and pollen flow (Provan 
et al., 2001). Therefore, they are useful for analysis of popula- 
tion genetics, genetic diversity, paternity analysis, and germ- 
plasm resource identification (Provan et al., 2001). In this study, 
we developed 12 cpSSR markers for V. unguiculata and evalu- 
ated their transferability to a related legume species, P. vulgaris. 
These results will be helpful for the future exploration and germ- 
plasm conservation in both V. unguiculata and P. vulgaris, 
although chloroplast microsatellite diversity in P. vulgaris 
has been investigated (Angioi et al., 2009; Desiderio et al., 
2013). 



METHODS AND RESULTS 

The complete chloroplast genome sequence of V. unguiculata was down- 
loaded from GenBank (GenBank accession no. NC_018051). The cpSSR loci 
distributed throughout the V. unguiculata chloroplast genome were screened 
using SSRHunter 1.3 software (Li and Wan, 2005). SSRs were selected based 
on the length of the core repeat motif (>10 nucleotides), for example, five units 
of dinucleotide repeat motifs, four units of trinucleotide repeat motifs, or three 
units of tetranucleotide repeat motifs. Primer pairs were designed based on the 
flanking regions of each SSR locus using Primer3 (Rozen and Skaletsky, 2000). 
The parameters of each primer were set using the following criteria: (1) primer 
size of 20-24 nucleotides in length; (2) GC content of 40-60%; (3) annealing 
temperature between 50-60°C; and (4) expected amplicon size of 100-300 bp. 
In total, 15 cpSSR primer pairs of V. unguiculata were designed and synthe- 
sized (Sangon, Shanghai, China). Twelve of them showed polymorphic bands 
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Table 1. Characteristics of 12 polymorphic cpSSR markers developed in Vigna unguiculata. 



Repeat 

Locus motif Primer sequences (5'-3') T a (°C) 



in V. unguiculata accessions, two were monomorphic, and one primer pair gave 
no products. The 12 polymorphic markers were used in the following analysis. 

A total of 91 samples were used in this study, including 62 V. unguiculata 
accessions and 29 P. vulgaris accessions (Appendix 1). All the samples were col- 
lected from an agricultural field in Anshan (30.46°N, 113.94°E), Caidian Dis- 
trict, Wuhan City, and preserved in Hubei Province Engineering Research Center 
of Legume Plants, Wuhan, China. Tender young leaves of each sample were col- 
lected and stored at -80°C until use. Total DNA was extracted from all the 
samples using the cetyltrimethylammonium bromide (CTAB) method (Doyle and 
Doyle, 1987). The yield and purity of the DNA were measured using a spectro- 
photometer SP-1910UVPC (Shanghai, China) at an A260/A280-nm wavelength. 

Characteristics of cpSSR markers were examined in both V. unguiculata 
and P. vulgaris. The same PCR conditions were applied in the two species. The 
PCR amplifications were performed in a 20-jiL reaction mixture containing lx 
Taq buffer, 30 ng of genomic DNA, 1.5 mM MgCl 2 , 200 uM dNTPs, 0.5 \iM 
for each primer, and 0.5 U Taq polymerase (MBI Fermentas, Vilnius, Lithua- 
nia). The PCR conditions were as follows: an initial denaturation at 94°C for 
5 min; followed by 35 cycles of 30 s at 94°C, 30 s at the locus-specific annealing 
temperature (Table 1), and 40 s at 72°C; and a final extension at 72°C for 5 min. 
The PCR products were separated using 6% denaturing polyacrylamide gels 
(Acr : Bis = 19:1) and visualized with silver staining. Due to the nonrecombin- 
ing nature of the chloroplast genome, each pair of chloroplast microsatellite 
primers was considered as a "locus" at a cpSSR site. Length variants of chloro- 
plast microsatellites at each cpSSR site were treated as alleles. Alleles detected 
from polymorphic primer pairs were used to generate a chloroplast haplotype of 
each individual; multilocus haplotypes were obtained by combining alleles 
from all polymorphic loci. Based on the polymorphic cpSSR markers, the frag- 
ment size amplified from each locus was scored by referring to a 20-bp DNA 
ladder (TaKaRa Biotechnology Co., Dalian, China). The number of alleles (A) 
and unbiased haploid diversity index (h) per polymorphic locus were calculated 
using the software GenAlEx version 6.41 (Peakall and Smouse, 2006). To estimate 
the informativeness of each SSR marker, the polymorphism information content 
(PIC) was calculated using the formula described by Botstein et al. (1980). 

As shown in Table 2, the characteristics of the 12 polymorphic cpSSR loci 
are tested in 62 V. unguiculata samples. A ranged from two to four in V. 
unguiculata (average: 2.75), h ranged from 0.123 (VgcpSSR4) to 0.497 
(VgcpSSR5) (average: 0.240), and PIC ranged from 0.114 (VgcpSSR4) to 
0.369 (VgcpSSR5) (average: 0.211). 

The transferability of the 12 V. unguiculata cpSSR markers was assessed in 
a related species, P. vulgaris; parameters of genetic variation were evaluated in 



GenBank Size range in Size range in 

Position 21 Region accession no. V. unguiculata (bp) P. vulgaris (bp) 

)-220 190-196 

)-190 186-190 

i-305 265-355 

)-280 180-280 

)-202 190-202 

)-196 196 

1-156 160 

)-182 182-186 

1-186 168-170 

»-220 170-220 

t-210 210 

t-210 178-180 



29 P. vulgaris individuals (the P. vulgaris group) (Table 2). All of the 12 cpSSR 
markers were successfully amplified in the P. vulgaris group, and nine showed 
polymorphisms, with the exception of VgcpSSR7, VgcpSSR9, and VgcpSSR13, 
which were monomorphic markers. Therefore, it indicated that 75% of these 
markers can amplify polymorphic bands. In P. vulgaris, A ranged from one to 
two, with an average value of 1.75. For each cpSSR locus, h was between 0.000 
(VgcpSSR7, VgcpSSR9, and VgcpSSR13) and 0.529 (VgcpSSRIO and 
VgcpSSR14) (average: 0.312). The PIC value varied between 0.183 (VgSSR3) 
and 0.374 (VgcpSSR2, VgcpSSRIO, and VgcpSSRH) (average: 0.312). 



CONCLUSIONS 

Twelve polymorphic cpSSR markers were developed in V. un- 
guiculata and showed high transferability in P. vulgaris. Further 

Table 2. Characterization of the 12 cpSSR markers in V. unguiculata and 
their cross-species amplification in P. vulgaris. 



V. unguiculata group P. vulgaris group 



Locus 


A 


h 


PIC 


A 


h 


PIC 


VgcpSSRl 


3 


0.210 


0.196 


2 


0.323 


0.262 


VgcpSSR2 


3 


0.362 


0.303 


2 


0.516 


0.374 


VgcpSSR3 


2 


0.153 


0.139 


2 


0.212 


0.183 


VgcpSSR4 


2 


0.123 


0.114 


2 


0.380 


0.298 


VgcpSSR5 


2 


0.497 


0.369 


2 


0.467 


0.332 


VgcpSSR7 


2 


0.125 


0.116 


1 


0.000 




VgcpSSR9 


2 


0.151 


0.138 


1 


0.000 




VgcpSSRIO 


4 


0.256 


0.237 


2 


0.529 


0.374 


VgcpSSRl 1 


3 


0.202 


0.185 


2 


0.441 


0.329 


VgcpSSR12 


4 


0.270 


0.255 


2 


0.349 


0.280 


VgcpSSR13 


3 


0.154 


0.146 


1 


0.000 




VgcpSSR14 


3 


0.383 


0.328 


2 


0.529 


0.374 


Average 


2.75 


0.240 


0.211 


1.75 


0.312 


0.312 



Note: A = number of alleles for each locus; h = unbiased haploid diversity; 
PIC = polymorphism information content. 



gcpSSRl (TA) 5 
gcpSSR2 (AAT) 5 



gcpSSR3 
gcpSSR4 
gcpSSR5 
gcpSSR7 
gcpSSR9 



(TA) 12 
(AT) 5 
(TA) 5 
(TA) 6 
(TA) 5 



gcpSSRIO (AT) 5 

gcpSSRl 1 (AT) 6 

gcpSSR12 (AT) 6 

gcpSSRl 3 (TA) 5 

gcpSSR14 (AT) 5 



GGTGGATGTTTATACCCAATCG 

TCTTTC T GC G AT AC AAAC AAG AA 

TTTTCTATGTATGGCGCAACC 

CGGGGATAAAGCTGCCTATT 

AAAC C AC T C GAAT AT T AT GGAAA 

CCAGTTCAAATCTGGTTCCTG 

GAAAAGAACAAGCAAATCCACA 

TGATCCTTACGATGCTTCCTTT 

AGCCCACTTTTCCGTAGGTT 

CTTTTCCTTGCCATAATGGTT 

TCAACCATTTCCCAACACCT 

CATCGAGTTCATGGATTTGC 

T G A A A T T T G A A A A A C G G G G T A 

AAGCGATACGGATAGATTCCT 

GGGCTCATTGGCTGTAGAAA 

CCATCTCTCCCCAATTGAAA 

TTTGAGAAGGTTCAATTGTTCG 

T CGGAC T C T AGGAAAGGAC AA 

GGCCATTTATCCCACTTTCC 

CCAGTCTCTACTGGGGGTTA 

TATTGGTTTTGCACCAATCG 

ACCAGGGTGTATGTGCGACT 

TGGATCATAATCCTTGAACATCA 

T GC GAAAAC AAAGAT AAGAAAT C A 



60 trnK-rbcL IGS LSC KF662476 190 

60 rbcL-atpBIGS LSC KF662477 180 

57 ndhJ-trnFIGS LSC KF662478 185 
60 ycf3 exon LSC KF662479 180 

58 P saB-rpsl4lGS LSC KF662480 190 

59 psbD-trnT IGS LSC KF662481 136 
57 trnR-trnSIGS LSC KF662482 144 
59 trnR-trnSIGS LSC KF662483 150 

59 petA-psbJIGS LSC KF662484 168 
56 psbJ-psbL-psbFIGS LSC KF662485 162 

60 rpl20-rpsl2 IGS LSC KF662486 162 
59 psaC-ndhEIGS SSC KF662487 162 



Note: IGS = intergenic 
a Position of each SSR 



spacer; LSC = long single-copy region; SSC = short single-copy region; T a = annealing temperature, 
in chloroplast complete genome of Vigna unguiculata (GenBank accession number: NC_018051). 
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analyses indicated that the cpSSR markers of V. unguiculata 
could reveal a relatively high level of genetic diversity in both 
V. unguiculata and P. vulgaris germplasm. These markers can 
be used to investigate genetic diversity and evolution in V. un- 
guiculata and P. vulgaris. 
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Appendix 1. Voucher information for legume species used for the cpSSR polymorphism study. All vouchers are deposited at the Hubei Province Engineering 
Research Center of Legume Plants, Wuhan, China. 



Species 


Voucher accession no. 


Country of origin/source 


Species 


Voucher accession no. 


Country of origin/source 


Phaseolus vulgaris 


B48 


China 




B4 


China 




C-l 


China 




B5 


China 




C-2 


China 




B6 


China 




C-3 


China 




B7 


China 




C-4 


China 




B8 


China 




C-6 


China 




B9 


China 




C-7 


China 




BIO 


China 




C-8 


China 




Bll 


China 




C-ll 


China 




B12 


China 




C-12 


China 




B15 


China 




(13*20)-2 


China 




B16 


China 




(13*20)-5 


China 




B17 


China 




(13*20)-10 


China 




B18 


China 




(13*20)-7 


China 




B20 


China 




(13*20)-1 


China 




A80 


China 




(13*20)-9 


China 




A89 


China 




(13*20)-4 


China 




A93 


China 




(1*7)-1 


China 




A96 


China 




(l*7)-2 


China 




A98 


China 




(1*7)-/ 


China 




A104 


China 




(l*7)-9 


China 




A105 


China 




(1*7)-10 


China 




A115 


China 




(l*7)-3 


China 




A125 


China 




(3*10)-4 


China 




A136 


China 




(3*10)-5 


China 




A143 


China 




(3*10)-6 


China 




Al 


China 




(3*10)-7 


China 




A8 


China 




(3*10)-8 


China 




A27 


China 




(3*10)-9 


China 




A33 


China 


Vigna unguiculata 


B28 


China 




A58 


China 




B30 


China 




A70 


China 




B32 


United States 




A156 


China 




B34 


China 




A162 


China 




B35 


China 




A168 


China 




B36 


China 




A171 


China 




B37 


China 




A176 


China 




B39 


China 




A181 


China 




B42 


China 




A182 


China 




J2 


United States 




A185 


China 




J3 


United States 




A189 


China 




J5 


United States 




A192 


China 




J7 


United States 




A194 


China 




J9 


Africa 




A71 


China 




Jll 


Mexico 




A82 


China 




Jl3 


Germany 




A155 


China 




B3 


Japan 
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